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Abstract 

This paper presents an analysis of the motion of an neutral impurity species 
in a nanometer scale 4 He cluster, extending a previous study of the dynamics 
of an ionic impurity. It is shown that for realistic neutral impurity-He po- 
tentials, such as those of SF 6 and OCS, the impurity is kept well away of the 
the surface of the cluster by long range induction and dispersion interactions 
with He, but that a large number of 'particle in a box' center of mass states 
are thermally populated. It is explicitly demonstrated how to calculate the 
spectrum that arises from the coupling of the impurity rotation and the center 
of mass motion, and it is found that this is a potentially significant source 
of inhomogeneous broadening in vibration-rotation spectra of anisotropic im- 
purities. Another source of inhomogeneous broadening is the hydrodynamic 
coupling of the rotation of the impurity with the center of mass velocity. A 
quantum hamiltonian to describe this effect is derived from the classical hy- 
drodynamic kinetic energy of an ellipsoid. Simple analytic expressions are 
derived for the resulting spectral line shape for an impurity in bulk He, and 
the relevant matrix elements derived to allow fully quantum calculations of 
the coupling of the center of mass motion and rotation for an impurity con- 
fined in a spherical He cluster. Lastly, the hydrodynamic contribution to the 
impurity effective moment of inertia is evaluated and found to produce only 
a minor fractional increase. 

I. INTRODUCTION 

The last few years have seen dramatic advances in the spectroscopy of atoms and 
molecules attached to large He clusters [1]]. These clusters provide a unique environment 
for a spectroscopy which combines many of the attractive features of both high resolution 
gas phase spectroscopy and traditional matrix spectroscopy [0. These include the abil- 
ity to obtain rotationally resolved spectra of even very large molecules such as SF 6 |§ and 
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(CH 3 )3SiCCH || (though with effective rotational constants only 1/3 - 1/5 as large as those 
of the gas phase molecule), and the ability to form and stabilize extremely fragile species ||, 
including high spin states of molecules. Particularly interesting is the recent demonstration 
by Grebenev, Toennies, and Vilesov |J that the free rotation is a direct consequence of 
the Boson character of 4 He, and can be viewed as a microscopic version of the famous An- 
dronikashvili experiment [0] that was used to measure the superfluid fraction in bulk liquid 
He. 

Despite rapid progress, many fundamental questions remain about spectroscopy in this 
environment. One such important question is the origin of the linewidth observed in ro- 
vibrational spectra. These linewidths vary with species, from 150 MHz reported for the R(0) 
line of the OCS z/3 fundamental band ||, to 1.5 cm -1 reported for H2O ||. Very recently, 
pure microwave spectra have been observed for impurity molecules (CH3CN and HCCCN) 



which have linewidths on the order of 1 GHz ||10|| . Since this is comparable to linewidths 
observed in ro-vibrational transitions, it is unlikely that vibrational dephasing, which is the 
dominant line broadening mechanism in the spectra of impurities in classical liquids, also 
plays a dominant role in He cluster spectroscopy. Likewise, vibrational population relax- 
ation, which has been invoked as the source of line broadening, obviously cannot contribute 
to the broadening in the case of microwave excitation. 

Given the fact that the He clusters are uniquely fluid, even down to zero temperature, it 
is natural to suppose that the impurity spectra will not display inhomogeneous effects. Vari- 
ations in local binding sites tend to dominate the linewidth of spectra in low temperature 
crystal or glass hosts. In a liquid environment, local changes in solvation lead to dynamic 
fluctuations in the spectral intervals, and thus can lead to dephasing but not static inhomo- 
geneous effects. Since the time scale for solvation fluctuations, due to the large zero point 
kinetic energy of He atoms in the bulk, is expected to be much faster than the dephasing 
times observed in most ro-vibrational spectra, one can expect the effects of the fluctuating 
solvation to be strongly motionally averaged. More fundamentally, it is far from clear that 
we expect pure vibrational dephasing at all given the low temperature of He clusters. Tak- 
ing the model of a crystalline solid, one expects dephasing to arise from interactions with 
thermal excitations. Since, as will be discussed below, the only excitations thermally excited 
are surface waves, and these are expected to interact extremely weakly with impurities, it 
is difficult to see how dephasing can be important. 

The fact that the existing experiments produce a broad distribution of cluster sizes 
suggests an inhomogeneous contribution to the line broadening from a size dependent shift 
in the absorption. However, it is likely that this effect is also not dominant in most spectra. 
First, the over all shift of ro-vibrational line from the gas phase to He clusters is fairly small, 
typically at most a few cm -1 0], so the spread in shifts as one approaches the bulk limit 
is expected to be small. Below, an estimate of the expected shift of the SFq z/3 vibrational 
transition due to the finite size of the He cluster will be given. Further, the experimentalist 
has some control over the mean of the cluster size distribution. In previous work on SF 6 , it 
was found that the spectral shift and temperature was essentially constant for He clusters 
of > 10 3 He atoms [|. 

An additional source of inhomogeneous environment needs to be considered, the distri- 
bution of the position and orientation of the impurity relative to the surface of the cluster. 
That this could be relevant is suggested by the following simple calculation. Consider an 
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impurity with an effective mass (including hydro dynamic contributions) of 50 atomic mass 
units (u) in equilibrium with a He cluster, which is known to be at 0.38 K ||. This molecule 
will have an root mean squared (RMS) velocity of 14 ms" 1 . In a cluster of diameter 10 nm, 
which corresponds to ~ 11,000 He atoms, it will take m 0.7ns for the impurity to cross 
the cluster. If the scattering of the impurity from the surface of the cluster leads to re- 
orientational dephasing, the resulting line broadening will be ~ 500 MHz, comparable to 
what has been observed in ro-vibrational transitions |]J. An additional source of dephasing 
is the hydrodynamic coupling of the molecular orientation to the molecular velocity, which 
will be described in more detail below. 

The simplest treatment for the motion of the impurity molecule is as a 'particle in a box'. 
This was considered in an earlier paper by Toennies and Vilesov [Tl||. That treatment was 
unsatisfactory on several grounds however. First, as we will demonstrate below, the true 
potential is not really 'flat' inside the box relative to the low value of the thermal energy. 
Second, and perhaps more importantly, they did not attempt to calculate the intensity of 
transitions in this model. Since the effective potential is almost the same in the upper and 
lower vibrational states, there are small Frank-Condon factors for transitions that change 
the center of mass quantum numbers of the impurity, and as such, changes in energy of the 
particle in a box levels are not directly observable in the spectrum of a neutral impurity. 

This paper will present a model that allows calculation of the effective potential for the 
motion of an impurity that is free to move about in a spherical liquid He droplet of radius 
R. A thorough analysis of the liquid drop model for pure He clusters was provided by Brink 
and Stringari |l2j]. The present calculation will be based upon long range interactions and 
treat the He as a continuum incompressible liquid of constant number density, p, which we 
will take as the known value of bulk liquid He. While such a model cannot be expected 
to correctly describe the interaction of the impurity with the first few solvation shells, this 
solvation structure is expected to be effectively constant as long as the impurity does not 
too closely approach the surface of the cluster. For realistic values of the impurity-He long 
range interaction, the energetics are effective at keeping the impurity from this difficult 
region. Further, since the impurity is kept away from the edge of the cluster, the fact that 
the He cluster has a diffuse (~ 7 A) interface with vacuum |13| will not be expected to 
significantly effect the impurity potential. Further, we can, with good approximation, keep 
only the lowest order terms in the expansion of the long range He-impurity interaction in 
inverse powers of distance. 

The paper is organized as follows. Section |IT] will consider a neutral, isotropic impurity, 
with SF 6 used as an example. Sections |I| and |IV| will consider the case of a symmetric 
top or linear molecule impurity, with HCN used as an example. Section [III] will derive 
the effective potential and apply a classical model for calculation of the resulting spectra. 
Section |IV| will derive the matrix elements needed for a fully quantum treatment of the 
interaction of the center of mass motion and the rotations of the impurity for a linear 
molecule, leading to a calculation of the expected spectral structure for the case of the HCN 
R(0) rotational line. Section [V] will consider the effect of the interaction of an impurity with 
the internal modes of the cluster, deriving the expected shifts produced by the interaction 
with the ripplons, which are quantized surface capillary waves. These interactions are found 
to be completely negligible contributions to the linewidth of transitions. Section [VI| will 
consider Hydrodynamic effects and derive the effective hamiltonian and matrix elements for 
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the hydrodynamic coupling of the rotation and center of mass momentum of the impurity. 
When this model is applied to the HCN R(0) line, poor agreement is found with experiment. 
However, it will be shown in Section |V11| that the combination of the anisotropic potential 
terms and hydrodynamic coupling leads to a prediction for the lineshape of the R(0) line in 
the OCS ro-vibrational spectrum that is in excellent agreement with that observed in Helium 
clusters and reported by Toennies Group Hl4 |. Section [VII1| will examine the hydrodynamic 
contribution to the reduction in rotational constant for an impurity, treated as an ellipsoid. 
It is found that the hydrodynamic kinetic energy only makes a minor contribution to the 
observed increase in moment of inertia of impurity molecules. Lastly, section [TX] will provide 
an overview and summary of this work. 



II. ISOTROPIC NEUTRAL IMPURITY, SUCH AS SF 6 , IN A HELIUM CLUSTER 



In a previous paper |]T3j| , the potential of an ionic impurity in a helium cluster was derived 
from the ion-induced dipole interaction of the impurity in Helium. The present paper will 
consider the case of a neutral impurity. As will be shown, this potential is dominated by the 
leading power in the long range interaction between the impurity and atomic helium. In the 
present section, the case of an impurity that has an isotropic long range interaction with He 
will be considered. This will be the case for impurity molecules of tetrahedral or octahedral 
symmetry as well as for atomic impurities. In the next section, the case of an anisotropic 
long range interaction, as expected for a linear molecule impurity, will be treated. 

As in our previous paper, we will take as a reference of energy not an isolated impurity 
atom or molecule, but the impurity in bulk He with constant He number density, p. In 
this way, the difference between an impurity inside the cluster and our reference state is the 
absence in the later case of He outside the cluster. We will also implicitly assume that the 
impurity is far enough away from the surface of the He cluster that we treat the interaction 
of the impurity with this 'missing' He as pairwise additive and we can use only the long 
range terms of the potential. 

Let the interaction energy of the impurity with the He be a van der Waals attraction at 
long range, with a long range form of : 

Where Cq and C§ are the long range dispersion interaction constants and r is the distance 
from the center of mass of the impurity to the He atom. For SF 6 , the lowest order anisotropic 
term appears at the r~ 10 level. 

Consider the impurity to be displaced from the center of a He cluster of radius R by a 
distance a. The z axis is defined as parallel to displacement vector of the impurity from 
the center of the cluster, which provides the origin for the coordinate system. Let r(6) be 
the distance from the impurity to the droplet surface at a polar angle 9 measured from the 
impurity. Basic trigonometry gives: 

R 2 = ( a + r cos(#)) 2 + r 2 sin(#) 2 = a 2 + r 2 + 2ra cos(#) (2) 

from which we can derive: 
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r (0) = JE 2 - a 2 sm(6) 2 - acos(6) 



(3) 



We can 'sum' up the missing interaction for all the He atoms absent (compared to the bulk 
phase) due to the finite extent of the cluster, using the interaction potential given above in 
Eq. 0: 
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We will consistently use y = a/R in this paper and refer to y as the normalized displacement 
of the impurity. We start with label F 2 to avoid confusion with the function F\ defined in 
the paper on ionic impurities. 

In atomic units, the C& and Cg coefficients are typically of the same magnitude ||16|| , 
and thus we expect V^ to be smaller than V® by a factor on the order of a^/R 2 ~ 10~ 4 
(ao is the Bohr radius), and thus we will neglect it. Using the Cq coefficient for SF 6 -He 
equal to 35 Eha[5 = 3.35 ■ 10~ 78 J m 6 |L7] (E h is the Hartree, the atomic unit of energy), p = 



0.022 atoms A (the density of liquid He at low pressure and temperature), and R = 3 nm 
(which corresponds to the size of a cluster containing w 2500 He atoms), we calculate V 2 ° = 
0.572 cm -1 • he ~ 2kbT c . kb is Boltzmann's constant and T c = 0.38 K is the temperature of 
nanometer scale He clusters. If we assume that the effective mass for translational motion 
of the impurity is M eS = 186 u (the SF 6 plus 10 He atoms to include hydrodynamic effects 
and some 'sticking' of solvated He atoms), we get a harmonic frequency for vibration of 
the impurity about the center of the He cluster of 790 MHz. Figures |X| and |IX| show the 
energy levels of this potential below 1cm -1 w Akt,T c , calculated using the Numerov-Cooley 
method |1~8| . Inspection of the Energy as a function of n (radial) and L (angular momentum) 



quantum numbers shows that the thermally well populated energy levels have a spectrum 
much more like that of a particle in a three dimensional harmonic well than a particle in a 
spherical 'box'. Figure |IX| contains a plot of the probability density for finding the impurity 
a distance a from the from the center of a R = 3nm cluster, calculated both from the 
set of eigenf unctions and from the classical Boltzmann distribution. We see an excellent 
agreement of the two, confirming the classical character of this motion even at T c . We see 
that the probability density is largely localized in the region with a less than about half 
the cluster radius. The RMS value for the displacement of the impurity from the cluster 
center is 0.96 nm. The highly classical nature of the thermal motion can also be seen by 
comparing the thermodynamic functions calculated from the exact quantum energy levels 
and for classical motion with the same potential. This comparison is shown in table 1. Also 
included in this table is the mean squared orbital angular momentum (calculated in units of 
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h 2 ), calculated as the thermally averaged value of L(L + 1), where L is the orbital angular 
momentum quantum number of the impurity motion around the center of the cluster, which 
is compared to the classical thermal mean value of < L 2 >= 2M eS k b T c < a 2 > /h 2 . The 
table also contains the thermodynamic averages for a cluster of R = 5nm as well as the 
corresponding properties of the ripplons (to be discussed below). It can be seen that the 
average energy and heat capacity are intermediate between those of a true particle in a box 
and a particle in a Harmonic well. The partition function for translation grows faster than 
the volume of the cluster, a factor of 15.3 compared to 4.6 when the radius of the cluster is 
increased from 3 to 5nm. This reflects the fact that in a larger cluster, the impurity gets 
to thermally sample a larger fraction of the cluster volume, as demonstrated by the RMS 
thermal displacement, \fa?. 



III. SYMMETRIC TOP OR LINEAR MOLECULE IMPURITY IN THE HELIUM 

CLUSTER 

We will now assume we have an impurity with an anisotropic interaction with He. If the 
impurity is in the center of the cluster, it will experience an isotropic net potential, due to 
the symmetry of the distribution of He atoms. However, as the impurity moves away from 
the center, it will experience an anisotropic potential. This anisotropy will cause a coupling 
between the center of mass motion of the impurity and its rotational motion, and can be 
expected to result in orientational dephasing. It is of interest to point out that the thermal 
mean value of the orbital angular momentum (L) of the impurity in the cluster « 20 H is 
much greater than the thermally populated value of rotational angular momentum of the 
impurity, (j), and as a consequence the procession of j about the total angular momentum 
J would be expected to completely dephase the average dipole moment, and lead to line 
broadening. 

The long range interaction of a neutral, symmetric top impurity with He can be written 
as the sum of two parts, one due to induction and the other dispersion. The former is 
dominated by the dipole induced dipole interaction for an impurity with a permanent dipole 



moment. We can write the potential in the following form [19 



V(r,4>) = ~- ^P 2 (cos(0)) - ^(008(0)) - ^P 3 (cos(0)) (9) 

where P n are the Legendre polynomials, and <ft is the angle between the vector from the 
center of mass of the impurity to the He, and the symmetry axis of the impurity. The 
induction contributions to these coefficients can be calculated |19| from the dipole moment, 
fj,, and quadrupole moment, 0, of the impurity and the polarizability, a, of He: 

2 

^<ind ^ind _ ^ ( T n\ 

£(ind _ ^ ^i n d _ 12 Q /i. Q , . 

73 ~ 3 71 5 (47re ) 2 1 ' 

The dispersion contributions to the long range interactions can be written as integrals over 
the imaginary frequency dependence of the polarizability of He and the polarizability and 
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dipole-quadrupole polarizabilities of the impurity [20] . Polarizablility is higher parallel than 
perpendicular to bonds. As a result, for linear molecules and most prolate symmetric tops, 
Cq2 > 0, which means that for such impurities, the interaction is lower in energy when the 
He-impurity displacement is parallel rather than perpendicular to the impurity symmetry 
axis. 

As an example, we will consider the case of HCN as the impurity. Using the parameters 
a(He) = (4 tc e ) -1.383 a|j, /x(HCN) = 1.174 e a = 2.98 Debye, and 6(HCN) 
we get that the induction contribution to the long range coefficients are: C\ 
1.14 eVA 6 and C^f 
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have been calculated by Atkins and Hutson |L9| to be Cgo Sp = 7.81 eVA 
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To find the effective potential of the impurity in the cluster, we will once again integrate 
the above effective long range interaction (Eq. |9|) over the He atoms 'missing' from the 
cluster. As before, we define the z axis as along the displacement of the molecule from the 
center of the cluster. We also define x as t ne angle between symmetry axis of the impurity 
and this z axis. In terms of these coordinates, we have the angle, 0, between this impurity 
symmetry axis and the vector from the impurity to a point in space as: 



cos(0(r')) = sin(#) cos(</>) sin(x) + cos(#) cos(x) 

where 6 and <p are the spherical coordinates of the point r' defined relative to the 
of the impurity. The results of the integration give the following results for the 
potential for the impurity: 
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F 2 (y) and V 2 ° are unchanged from the previous definitions. Despite the rather 
appearance of F±(y) and F e (y) near y w 0, they are regular: 



F,(y)* 3 -y 2 + ^ + ^y 6 ... 
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F 5 (y)^y + V + 10y 5 ... (21) 
F 6 (y)^^y 3 + ^-y 5 ... (22) 

If we consider the case of HCN in a He cluster of R = 3 nm, and use the parameters 
given above, we get for the effective potential: 

AE(a, R, x )^ = ( 0.2519F 2 (y) + 0.0045F 5 (y) cos(x)) 



he 



+0.0619F 4 (y)P 2 (cos(x)) + 0.0031F 6 ( 2 /)P 3 (cos(x)) ) cm- 1 (23) 



We see that the two C-j terms (F 5 (y) and F e (y)) can be neglected in first approximation 
compared to the isotropic and leading anisotropic term (F 2 (y) and F±(y)) that derive from 
the Cq terms. These coefficients are comparable to the line widths observed in ro-vibrational 
spectra of molecules inside of He clusters. The V 2 ° isotropic term will confine the impurity 
near the center of the cluster. As discussed above, this term determines the effective potential 
for the translational motion of the the impurity, while the weaker anisotropic term will couple 
the translational motion to the rotation of the impurity molecule. Figure IIX| shows a plot 



of these two terms as a function of a for the case considered here, along with a plot of the 
thermal probability density distribution as calculated from the isotropic potential, V^i^y). 
It can be seen that the anisotropic term, V^F^y), is much smaller than the isotropic term 
over the full range of displacements that are thermally well populated. 

Since Ceo > 0, it is lower in energy for the HCN axis to be perpendicular to the displace- 
ment from the center of the cluster (0 = tt/2). This is easily rationalized if one considers 
the case of HCN below a flat He surface. When cf) = tt/2, the nearest 'missing Helium' is 
perpendicular to the HCN molecular axis, and this cost less energy than when = 0. If the 
interaction energy is much weaker than the seperation of rotational levels of the impurity, 
which will be the case for HCN for all the thermally well populated configurations, to first 
order we can consider the effective potential for translational motion to be the average of 
the effective potential over the rotational function of the impurity. For a linear molecule, the 
odd Legendre terms average to zero, while the P 2 term will split the rrij degeneracy of each 
j level, such that the center of gravity is conserved. If we consider the impurity center of 
mass displacement as static, then the natural axis of quantization of the angular momentum 
is parallel to this displacement, with an energy of a j, rrij rotational state equal to: 

AE(a, J, M) = V 2 °F 2 (a/R) + V°F,(a/R) [ ^j - l)(2j + 3) ) ^ 

which implies that rrij = states have the highest energy and rrij = j the lowest of each 
manifold of fixed j states. 

A proper quantum mechanical treatment of the coupling of the rotation and translational 
motion of an anisotropic impurity will be given in the next section. This section will close 
by demonstrating a simple 'classical' application of the above results to estimate a line 
broadening mechanism in the case of SFg. The ground state of SFq has octahedral symmetry, 
and thus has only an isotropic long range interaction, as discussed above. However, the IR 
allowed vibrational excitations are to triply degenerate levels. The u 3 mode observed in 
the IR spectrum of SF 6 in He clusters has a particularly large vibrational transition dipole 
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moment of fi 3 = 0.388 Debye for the v 3 = — > 1 transition. As the SF 6 is displaced in the v 3 
mode, it will have a dipole moment proportional to normal coordinate Q 3 . This dipole result 
in a dipole-induced dipole interaction with the He, which modifies the force constant for 
this normal mode The modification of the v 3 wavenumber (relative once again to SF 6 
in bulk liquid He), as a function of the position of the SFg in the cluster and the orientation 
of the triply degenerate z/3 mode are given by Equations [H^ and [L3|, with the derivative of 
the dipole with respect to the dimensionless normal coordinate (which is just \[2 larger than 
the fundamental transition dipole moment) replacing the permanent dipole moment. If we 
treat these coordinates classically, we can calculate the spectrum of the v 3 mode of SF 6 in 
the cluster, relative to SF 6 in bulk liquid He from the effective potential. Thus, we can write 
z/ 3 (a,x) as 



v 3 {a, R, X ) = ^s(oo) + z/° [F 2 (a/R) + F A (a/R)P, 
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(25) 
(26) 



^3(00) is the frequency of z/3 in the bulk He limit and x is the angle between the diplacement 
of the SF 6 and the direction of polarization of the v 3 mode. Note that this expression would 
be exact if a did not change in time, i. e. we have a static inhomogeneous distribution 
of impurity positions, as in a matrix without diffusion. Dynamic motion of a and x will 
tend to motionally average the instantaneous frequency u 3 (a(t), x(t)) an d lead to a narrower 
spectrum than predicted by this approximation, which is sometimes known as the 'statis- 
tical method' and was first introduced by Kuhn and London [23| to model pressure 
broadening. 

For the case of the u 3 mode of SF 6 in a R = 3 nm He cluster, v° = 16 MHz. Since the 
isotropic term, F 2 (y) dominates over the anisotropic term, F^y), over the range of impurity 
displacements that are thermally populated, we will keep only the isotropic term. We can 
write the spectral density, S(v) as 
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V2F 2 (a/R)/k b T c )da. Figure [IX| shows the calculated spectral 



where Z c = J Ana exp 
lineshape for v 3 of SF 6 in this model. It can be observed that the residual shift of the 
mode (compared to the limit of SF 6 in bulk He) and the FWHM are both estimated to 
be about 30 MHz. This can be compared to the linewidth of ~ 300 MHz FWHM observed 
by Hartmann et al [Q. for the v 3 mode of SF 6 in 4 He clusters of mean size 4000 atoms 
(R = 3.5 nm). It is clear that the calculated spectral broadening is more than an order of 
magnitude below that observed in the laboratory, which argues that this effect contributes 
at most a minor fraction of the total dephasing rate. 

In their earlier paper that considered 'particle in a box' states of an impurity in He 
clusters, Toennies and Vilesov Oj proposed that these levels could be directly observed as 
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satellite lines on a vibrational band of an impurity. They suggested that such particle in a 
box state changing transitions could gain intensity due to differences in the effective potential 
in the upper and low vibrational states, but reported no estimate of the expected intensities 
of such transitions. Using the SF 6 example, the author has calculated Frank-Condon factors 
for transitions that change the quantum numbers for the center of mass motion along with 
a vibrational transition. The largest such Frank-Condon factor was below ~ 10 -4 . This 
is not unexpected since the excited state has an effective potential that is stronger by only 
~ 0.1%. Based upon this analysis, it appears unlikely that such center of mass motion 
changing transitions will be directly observed unless the change in potential with vibrational 
coordinate is much larger. It is worth pointing out that the spectral structure that Toennies 
and Vilesov suggested may be due to such transitions was latter assigned by them to the 
rotational structure of the vibrational transition inside the He cluster H. 



IV. INTERACTION OF IMPURITY TRANSLATIONAL AND ROTATIONAL 

MOTIONS 



Because the interaction between the impurity and the He cluster depends upon both the 
position of the impurity in the cluster and the orientational angles of the impurity, there is 
a coupling that can lead to relaxation of the orientational order, and thus line broadening in 
the impurity spectra. The purpose of this section is to elaborate how this can be calculated. 
While the physical origin of the effect is quite different, this problem has some mathematical 
similarity to the coupling of rotation and translation that has been used to model the IR 
spectra of light rotors in solid noble-gas matrices |23-|27| 



We will assume our impurity is free to move in the cluster with an effective mass M c g. 
We will further assume it to be a linear molecule with an effective rotational constant _B e fr 
in Helium. We will model the system with the hamiltonian: 

H = y2 + B ^ ' i 2 + V " F M + v2Hv) ■ ^(cos(M) (so) 

In this equation, V 2 operates on the center of mass coordinates of the impurity, j is the 
angular momentum operator for rotation of the impurity and 8u is the angle between the 
displacement of the impurity from the center of the cluster and the direction of the linear 
axis of the impurity. The results are easily extended to include other Legendre terms, such 
as those that arise from the r~ 7 terms in the long range potential of an isolated impurity 
and He atom. 

If we keep all but the last term in our zero order hamiltonian, then rotational and 
translational motion are separable and we can write wavefunctions as: 

\n,j, mj ,L,M L >=^F J , mj (n 1 )F L ,M L (n 2 ) (31) 

where Xn,L(a) are the numerical solutions to the radial potential, U 2 °F2(?/), with orbital 
angular momentum quantum number, L. Vt\ = 9i,<pi are the orientational angles of the 
rotor, and fl 2 = #2, 02 the angles of the center of mass of the impurity, but referenced to 
some arbitrary lab fixed axis system. Introduction of the V®F i (y)P2 term will couple j and 
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.L, to form a total angular momentum J. We can write the basis function of the coupled 
representation as: 

\n,j,L,J,M >= (-iy- L V2J + T £ ( ^ ^ _ J M )\n,j,m 1 ,L,m 2 > (32) 



mi ,m2 



The matrix elements of the zero order hamiltonian are diagonal in this representation as 
well, and give by: 

<i,Z = <L + 5eff-jO' + l) ( 33 ) 

The E® L and Xn,i( a ) can be determined by numerical solution (using the Numerov-Cooley 
method) of the radial Schrodinger equation. 

The coupling matrix elements can be calculated using standard methods of angular 
momentum theory |2"E||2"5| |, including the Spherical Harmonic Addition Theorem to write 
P2(cos(8ij)) in terms of products of spherical harmonics of coordinates Qi and f2 2 : 

<n,j } L,J } M\H'\n' } j / ,L',J,M> = V° < n, L \F A (r)\n' , V > .{-l) J+ i +L ' x 

Wi + W + WL + m + 1) • ( i { o ) ' ( ) ■ { 2 i< f } < 34 ' 

The two 3 J symbols in this expression lead to the selection rules for coupling that j' = j, j±2 
and V = L, L ± 2. Explicit expressions for all the required 3 J and 6 J symbols can be found 
in the tables given by Edmonds P5| . Since for the cases considered at present, the coupling 
matrix elements are comparable to the separation of the center of mass motion states, but 
are much less than the seperations of the rotational energy levels of the impurity, the matrix 
elements that are off diagonal in j are neglected in numerical calculations presented below. 

The only change needed for other Legendre terms is to replace the "2" in the above equa- 
tions by the order of the Legendre term, and the radial integral will be over the appropriate 
function, F n (y). These terms will loosen the selection rule on Aj and AL for states that can 
be coupled. Parity, which is determined by (— 1) J , will remain a good quantum number. 

In order to calculate the dipole transition moments, we will treat the electric field as 
along the laboratory Z axis. This gives: 

(y ijmi |cos(0) \Y fjmi ) = V(2j + l)(2j' + l)(-ir ( o ) " ( 4n X o) (35) 

where 9 is the angle between the Z axis and the symmetry axis of the impurity molecule. 
Going to the coupled representation we have: 

{n,j,L,J,M\ooB[e)\n',j',L',J',M') = (_i)^"+^+ 1 5 n , n ,8 L , v 6 m ,m> x 

j f 1 \f J J' l\ j L f J' 



® + + ^ 2J + + ^ { J { -M M J • { 1 J j ) ^ 

We can expand the eigenf unctions in terms of the functions of the coupled representation 



as: 



\k,J } M>= \n,j,L } J } M><n,j,L,J\k,J> (37) 
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(k is a label to distinguish eigenstates with the same total angular momentum quantum 
numbers). We do not put the M label on the amplitudes since these are independent of 
M due to rotational invariance. We can define a line strength factor for the intensity of a 
transition between two levels k, J and k', J' as follows: 

S(k, J; k', J') = 3 I < k, J,M\ cos(9)\k', J',M' > | 2 (38) 

M,M' 

by exploiting the normalization of 3J symbols, we can reduce this to: 
S(k, J; k', J') = (2 J + 1) (2 J' + 1) x 



£ < k, J\n j3l L, J >< n,f, L, J'\k \ J' > ^{2j + \){2j> + 1) 

n,j,j',L 



-1 



\j+j'+L 



( J f 






f 




\0 


J) 




J 





(39) 



Figure [IX| shows a simulation of the j = — > 1 rotational transition of HCN in a R = 
3 nm He cluster. The calculated 'stick' spectrum has been convoluted with a Lorentzian 
lineshape with FWHM of 3 MHz. A spectral broadening of the transition of ~ 150 MHz can 
be seen in the figure. The residual 'fine' structure would be washed out when one does an 
average over the experimental cluster size distribution in a simulation of a real experiment, 
but the overall FWHM of the transition is not a strong function of cluster size. Also shown 
in the figure is the lineshape calculated from the semiclassical 'statistical' spectral model, 
where the spectrum is calculated as a Boltzman average over different 'static' center of mass 
positions for the HCN. The two 'bumps' correspond to J, M = 0, — >• 1,±1 (more intense 
low frequency side) and — ► 1,0 (broader and less intense high frequency bump). These 
individual transitions are 'washed out' in the full quantum spectrum because the motion of 
the impurity causes the natural quantization axis (the direction of the displacement from 
the center of the cluster) to average over all directions. The linewidth observed for the R(0) 
line in the first C-H overtone band of HCN is ~ 1. GHz, an order of magnitude larger. This 
suggests that another mechanism dominates the dephasing rate in the experiment. 



Recently, Nauta and Miller |30| have observed the R(0) transition in the C-H funda- 
mental band of HCN. They find that the lineshape consists of two features, with a width 
and seperation in qualitative agreement with the above 'statistical' model, i.e., as predicted 
without motional averaging. Furthermore, they have studied the spectrum as a function of 
external electric field, and have found that the two features correlate with the two high field 
m — ±1 and m = peaks expected for the R(0) transition of a polar linear molecule. 



V. INTERACTION OF IMPURITY WITH HE CLUSTER EXCITATIONS 

Up to now, we have considered the effect of the cluster potential on the rotational and 
translational degrees of freedom of an impurity. We will now consider the interaction of the 
impurity with the internal vibrational excitations of the He cluster. These can be separated 
into two types. The first are bulk excitations that consist of density fluctuations or phonons. 
The second are surface excitations that consist of capillary waves and are known as ripplons. 
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As discussed by Brink and Stringari [0, for a cluster of N atoms, the lowest phonon 



excitation, the symmetric breathing mode, has a frequency of 540 iV -1 / 3 GHz. Even for 
a cluster of N = 10 4 He atoms, this corresponds to 25 GHz which is much larger than 
kbT c /h = 8.0 GHz. Thus, only a few percent of such clusters have any quanta of excitation 
in the phonon modes. This implies that such modes cannot produce a significant level of 
inhomogeneous broadening in impurity spectra. It is possible to observe excitations of the 
phonon modes in an impurity absorption spectrum, but such excitations should appear as a 
well separated 'phonon side band' which peaks several cm -1 to the blue of the 'zero phonon 
line' which does not involve phonon excitations. Such phonon wings are commonly seen 
in the electronic spectra of impurities in He pl|j3^| , but have not yet been observed in 
impurity vibration-rotation spectra. This is presumably due to a small change in effective 
He solvation around impurities following vibrational excitation. As a result, almost all the 
transition intensity is concentrated in the zero phonon line. 

In contrast, the lowest frequency surface waves have much lower wavenumbers, and exci- 
tations are thermally excited even at the low temperature of He clusters. These excitations 
involve no change in He density and the only restoring force is due to the increase in surface 
area caused by the change in the surface shape. The theory of such surface waves is pre- 
sented in the classic text on Hydrodynamics by H. Lamb |B3], and discussed, in the context 
of He clusters, by Brink and Stringari [0 , who calculated the spectra and thermodynamic 
functions of these excitations as a function of cluster size. 

Due to excitation of ripplons, the surface of the He cluster develops oscillations. If we 
consider a point of the surface of the cluster at polar angles 9 s ,<p s (again, relative to some 
arbitrary lab fixed coordinate system), we can write the displacement of the surface as: 

r(0 B , V s)=R+ E S ( L ^ M s )Y LsMs (9 S} <^ s ) (40) 

L S ,M S 

In the above equation, S(L S , M s ) is the normal coordinate of the surface wave with angular 
momentum quantum numbers L s , M s . The theory assumes that S(L S , M s ) <C R. The sum 
is for L s > 2. The functions for L s = correspond to the breathing mode, and L s = 1 
corresponds to a translation of the entire cluster, not distortions of the cluster shape. If we 
assume that the surface tension of He is a = 380 /iJm -2 fl2]| , the frequency of the L s , M s 
surface mode is independent of M s and is given by: 



= = hf^^^m < 41) 



/i -( i -- 1 )( £ - +2 )sk JV " /2 



= 78 GHz y/L B (L e - 1)(L S + 2)N- 1/2 

For a cluster of size N = 10 4 , v{L s = 2) = 2.2 GHz. The clusters will have a maximum 
value of (L s ) max « (7r/2)iV 1//3 which is estimated by assuming that the smallest wavelength 
is ~ twice the mean interparticle distance. For a cluster of R = 3nm, this corresponds 
to (L s ) max = 26. Table 1 presents thermal average quantities for the ripplons in clusters 
of sizes R = 3 & 5nm, which can be compared with the thermal averages of the impurity 
center of mass motion presented. Note that the average energy and entropy in the ripplons 
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greatly exceeds that of the impurity and forms the dominant 'heat bath' of the cluster. 
However, the average thermal angular momentum in the impurity center of mass motion 
approximately equals the total in the ripplons. We can thus anticipate that in coming to 
equilibrium, the angular momentum constraints may play an important role. It is also worth 
noting that for a cluster of size R = 3 nm, the total canonical internal energy at T c is about 
equal to the energy required to evaporate a single He atom, pa 18.SkbT c fl^] . Thus, one 
can anticipate that smaller clusters will be slightly warmer than larger clusters, due to their 
low heat capacity. This is in contrast to the fact that we expect a slight decrease in He 
binding energy for smaller clusters, and thus higher vapor pressure, which would argue that 
smaller clusters should be colder. The decrease in binding energy, due to change in surface 
free energy, can be estimated as 2a /(pR) = 1.73 ■ cm -1 • nm/R, and is considerably larger 
than the increase in binding energy caused by attraction to a neutral impurity, which is 
pa 0.17cm _1 -nm 6 /i? 6 for the case of SF 6 . 

We can model the interaction between an impurity and a ripplon by considering the 
change in the sum of impurity - He interactions caused by the modulation in the surface of 
the cluster. Let 5r(9 s , <p s ) be the modulation in the cluster radius in the direction (9 S , ip s ) 
measured from the center of the cluster. (Note: this is a change from our previous treat- 
ments where we used polar angles defined from the impurity). As above, let y equal to the 
normalized distance of the impurity from the center of the cluster. To first order in 8r, we 
can write the change in He-impurity interactions as: 

C m p r 5r(Q s )dQ s C e2 p f 5r(Q s )P 2 (cos(9 

[V) = 1VJ [l + y 2_ 2 ycos(^)] 3 ~WJ [l + y2-2y C os(e is )] 3 1 J 

where Cqq and C$2 are the isotropic and anisotropic contributions of the Cg coefficients, Q s 
is the surface solid angle, 9{ s is the angle between the vector from the center of the cluster 
to the impurity and the vector from the center of the cluster to the surface, and 9 isr is the 
angle between the vector from the impurity to the point on the surface in direction 9 S , 4> s 
from the center of the cluster and the vector along the symmetry axis of the impurity. It 
is also useful to introduce the additional quantities as follows: r- ls , the distance from the 
impurity to the surface; 9- a , the angle between the vector from the center to the impurity 
and the vector along the axis of the impurity; and 9 sr , the angle between the vector from 
the center to the surface and the vector along the axis of the impurity. We can write these 
angles in terms of the three sets of polar angles for the impurity center of mass (9i, <pi), the 
surface (9 s ,ip s ), and the rotor (9 T ,(p T ) as follows: 

cos(6 l i S ) = cos(6' i ) • cos(6' s ) + sin(6 , i ) • sin(6 l s ) ■ cos(<£>i — ip s ) (43) 
cos(6'i r ) = cos(6'i) • cos(6 l r ) + sin(#j) • sin(0 r ) • cos(v?i — ip T ) (44) 

r is = ayjl + y 2 -2yco8(6 ]s ) (45) 

cos(0 isr ) = — [cos(6y - ycos(9 il ) ] (46) 

'is 

An exact treatment of this model would consider the simultaneous interaction of the 
surface waves, the center of mass motion of the impurity, and the rotation of the impurity, 
each of which has its own angular momenta that must be coupled. We will consider here 
instead a model for the interaction of the ripplons and impurity rotation alone. As such, 
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we will treat the impurity as fixed at some normalized distance, y, from the center of the 
cluster in a direction we will select for the z axis. In the end, we will predict the spectrum by 
averaging over the possible positions y, using a Boltzmann weighting. In this approximation 
we have cos(9 is ) = cos(# s ) and cos(# ir ) = cos(6> r ). In order to predict the j = — > 1 
rotational transition of the impurity, we will need to average the interaction over the lowest 
rotational functions of the rotor. This gives: 



(Y 00 (Q r )\P 2 (cos(e ist )\Y 00 (n r ))=0 



(Kio(fir) |P2(cos(0 isr )|Yio(O r )) = - 



3 cos 2 



Aycos(e s ) + 2y 2 -l 



1 + y 2 — 2ycos( 



(47) 
(48) 



(Y 1±1 (fi r ) |P 2 (cos(0 isr )|Y 1± i (fi r )> 



< Y 10 (Q T )\P 2 (cos(9 isi )\Y 10 {Q r ) > 



(49) 



We thus find that the j = rotor state does not couple to the ripplons, while the j = 1 
state does. 

In order to calculate matrix elements of the interaction, we need to determine the di- 
mensionless amplitude for the ripplon vibrations. Using the methods presented in Lamb's 
text we can determine the effective mass, fi(L 8 ), for each ripplon: 

M HeP R 3 



fi>{L 8 



(50) 



Note that the effective mass is proportional to the total mass of the cluster, not the mass 
of atoms near the surface. This is due to the fact that the surface waves have a velocity 
field that goes into the bulk of the cluster, with the velocity proportional to y Ls ~ x ■ Surface 
waves of higher L s values are increasingly confined to motion of atoms closer to the surface, 
and thus have a reduced effective mass. Using this relationship, we can write the surface 
displacement in terms of raising and lowering operators for the ripplons (a'(L B , M 8 ) and 
a(L s ,M s )) as follows: 



5r(6 s ,<p s ) = £ S(L s ,M s )Y LsMs (6 s ,<p s ) 

L S ,M S 

S(L B , M s ) = d s (L s ) ( a\L„ M s ) + a(L s , M s ) 
d s (L s ) - 



h 



2p(L s )u(L s 



LJi' 



2(L s -l)(L s + 2)aM UcP R 3 



1/4 



(51) 
(52) 
(53) 



For R = 3 nm, d s (2) = 0.37A. Summing the thermally averaged displacement of all ripplons 
up to (-Ls)max gives an RMS thermal motion of the surface of 1.34(1.50) A for clusters of 
R = 3(5) nm radii. 

We can thus write the interaction of ripplons with the rotation of the impurity molecule 
as follows: 



n 



ripplon, rotation 



U £ (a\L s ,M s )+a(L s ,M s 



x 



L s , Ms 



1/4 



(L„- !)(£„ + 2) 



Y LsMs (n s )p 2 (cos(e iST ))dn s 

[l + y 2 -2y C os(6 is )] 3 



(54) 
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C eo p 
R 4 



ti 2 



2 a M He pR 3 



1/4 



(55) 



For an HCN impurity in a 5 nm He cluster, Uq = 1.3 MHz • h 

Despite the fact that the integral diverges as y — *> 1, the coupling matrix elements 
diagonal in rotational quantum numbers are much less than the spacing between ripplon 
states for the range of y that includes most of the probability density of the impurity. 
Since the coupling is linear in the ripplon coordinates (which are harmonic oscillators), 
7^rip P ion,rotation has no first order contribution to the energy. This suggests that we can treat 
this coupling by second order perturbation theory to account for the virtual creation and 
annihilation of ripplons as an impurity rotates. A consequence of the linear ripplon coupling 
is that this second order correction is independent of the occupation of the ripplon modes! 
Such coupling shifts the equilibrium position of the ripplons from zero displacement but 
does not change their force constant. We find a shift in energy equal to: 



AE(y,j,mj) = 



II 2 



1/2 



Ls>Ms hu(L s )\(L s -l)(L s + 2) 
Y LsMs (n a )(Y j , m m i )\P 2 (cos( 



I 



n 2 



[l + J/ 2_22/cos(0 fa )] i 



(56) 



If we select the z along the vector from the center of the cluster to the impurity 

center of mass, then only the M s = terms contribute to the sum. 

The above expression was used to calculate the shift in the energy of the j = l,m,j — 
level as a function of y. The integrals over 6 S were performed numerically using Gauss- 
Legendre integration. The shift for the two j = l,rrij = ±1 levels will be just one quarter 
as large. Figure |TX| shows the calculated shift divided by Planck's constant as a function of 
a, along with a plot of the radial distribution function for the impurity. This was calculated 
from the isotropic, Ceo, interaction energy, as described above. It can be seen that the 
ripplon induced shift is completely negligible except for the extreme tail of the Boltzmann 
distribution. We can conclude that the potential interaction of the impurity rotation with 
the ripplons can be neglected compared to the position dependent anisotropic dispersion 
interaction. 



VI. HYDRODYNAMIC COUPLING OF IMPURITY MOTIONS 

While the motions of the impurity in the superfluid He are expected to be without 
friction, the medium does exert an influence by hydrodynamic effects. The reason is that as 
the impurity rotates or translates, He atoms must move out of the way to avoid the short 
range repulsion. For a spherical impurity, the effect of the He kinetic energy produced by the 
impurity translational motion can be modeled by an effective mass for the impurity that is 
increased by 0.5 times the mass of He displaced by the impurity. For an anisotropic impurity, 
there is in addition both a hydrodynamic contribution to the effective moment of inertia 
for rotation, and a coupling between rotation and center of mass motion. In this section a 
detailed calculation of these effects will be presented. This type of coupling was previously 
discussed (though developed in much less detail) in a paper by Elser and Platzman, who 



16 



proposed an anisotropic solvation structure around positive ions in bulk liquid He to explain 



an observed temperature dependence to the effective mass of ions [34 



The hydrodynamic coupling of translation and rotation does not occur for an impurity 

. Let us model the impurity molecule as a rigid, 
3. 2 A and c = £a = 3.7A. These parameters are 



with tetrahedral or octahedral symmetry |33 
symmetric ellipsoid of principal axes a = b = 
based upon the position of the inner wall in the He-HCN potential of Atkins and Hutson |]19" 



Since the dominant hydrodynamic effect is for parts of the fluid close to the moving body, the 
use of a continuum, incompressible fluid model is a questionable approximation. However, 
because of its simplicity, it is worth working out in detail the predictions of this model. 
Further, it may be the case that the form of the resulting hamiltonian may be more generally 
applicable than the simple hydrodynamic model used to derive it. 

The motion of an ellipsoid in an incompressible fluid without viscosity is worked out in 
the text by Milne- Thomson ||35|| . The hydrodynamic kinetic energy is given for an ellipsoid 
moving in any direction relative to its principal axes, v denotes the velocity of the ellipsoid 
relative to the fluid at large distance from the ellipsoid; A denotes a unit vector in the 
direction of the unique principal axis; M' = Mjj e p | n a 3 £ is the mass of He displaced by 
the ellipsoid (34.6 u for the HCN-He parameters given above), and M the mass of the 
ellipsoid. The classical kinetic energy of translational motion of the ellipsoid is: 



t=-(m+- 

2 V 3 



+ -M'(2 Tl + a x )j v-v + -M'(on 



-7iJ 



' A • v) v ■ v 

3 



where 



a i 



2 - a 



Of, 



Jo JT 



+ x)(e 2 + x) 3 / 2 



7i 
7o 



7o 



2-7o 



;i + x) 2 v / ft^ 



(57) 

(58) 
(59) 



Figure |TX| shows a>i, 71, and their difference as a function of £. We can see that the coupling 
between translational motion and rotation is zero in the case of a sphere (£ = 1), saturates 
for a very prolate ellipsoid (£ 1), and diverges for an extremely oblate ellipsoid (£ — ► 0), 
i.e. for a flat disk. For our model of HCN in He, £ = 1.16, and we have a± = 0.419 
and 71 = 0.544. The coupling of A and v reflects the fact that it takes less energy (at 
fixed velocity) to move the ellipsoid though the fluid with its long principal axis aligned 
with the flow, then with this axis perpendicular to the flow. This leads to coupling of 
the rotational motion to the translational motion. In the bulk, the rotational axis will 
precess about the velocity vector, and this would be expected to produce re-orient at ional 
dephasing of a molecular ro- vibrational absorption line, since there is a distribution of 
molecular translational velocity. 

In order to derive the quantum hamiltonian for the motion of such an ellipsoid, we need 
to re-express the kinetic energy in terms of the momentum of the ellipsoid, defined by: 

p = V V T =[M + M' 7l ] v + M' («! - Tl ) (A 

which can be inverted to write 



v) A 



(60) 



M + M'7i 



_ M'(ai-7i)(A-p) ^ 
P M + M'a ± 



(61) 
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Substituting into the above expression for the kinetic energy gives 



T = \V' V- -, r (A • p) (62) 

2(M + M'-f 1 y y 2(M + M'a 1 )(M + M / 7l ) V y J v ; 

Making the substitution p — > ihV we get the quantum hamiltonian: 

(A ■ V) 2 - ^V 2 | + B cS J 2 + V(f, A) (63) 



h 2 _ 9 „ f/r -\2 1, 



where 



3M + 2M'«i + M'7i V ; 

H 2 (M + Af'ttx) (M + M'71) 1 J 

For our model of HCN in He, C H A = -0.772 GHz A 2 , and M eff = 35.2 u. 

Before considering the case of an impurity in a cluster, it is instructive to examine the case 
of a linear molecule impurity in bulk He, where the momentum is a good quantum number. If 
we assume that we can treat the hydrodynamic coupling by first order perturbation theory, 
the shift in energy for a level j,rrij (with the quantization axis along the velocity), and 
translational kinetic energy -Etrans, is given by: 

E hydm (j, TTlj ) = A H £trans ( 2 j - l)(2j + 3) ^ ' 

v ; \ v 3M + 2M'ai + M / 7 1 y v ; 

For the case of HCN, A-^ = 0.036. The j = — > 1 transition will be shifted by +0.4An-E'trans 
and — 0.2y4n-Etrans for polarization parallel and perpendicular to the velocity vector of the 
impurity. The width of a j — ► j+1 transition will scale as Figure |IX| shows the calculated 
spectral lineshapes for the R(0) and R(2) transitions after averaging over the Maxwell- 
Boltzmann distribution for ^trans- The horizontal axis is in units of A-^k^Tch" 1 which equals 
283 MHz for the case of HCN. The spectra have a cusp in the center due to the \J E trans 
factor in the density of states. The higher frequency lobe arises from parallel transitions and 
the lower lobe, the perpendicular transitions (relative to the impurity momentum vector). 
Note that molecules with the projection of the rotational angular momentum perpendicular 
to the center of mass momentum are higher in energy than those with rotational angular 
momentum either parallel or antiparallel to the linear momentum. This is because in the 
constant momentum case, orientations with the linear molecule axis parallel to the linear 
momentum are higher in energy, due to the lower effective mass for that orientation. 

Explicit expressions for the matrix elements of the hydrodynamic term can be derived 
from the relationships from the matrix elements of A- V operating on the spherical harmonics, 
as given by Rose . After considerable algebra, the matrix elements diagonal in j in the 



coupled representation (where L is the center of mass motion orbital angular momentum 
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quantum number, j the quantum number of rotation of the linear impurity, and J the 
quantum number for the total angular momentum of the impurity), 



\n,L,j,J,M>= ]T (-l) L -i +M V2JTT 

mi,m2 

can be found. For V = L + 2 they are: 

,2 1 



L j J 
mi m 2 — M 



>L, mi (fi)^, m2 (^r) (68) 



<n',L + 2,j,J,M|(X-V) --V 2 |n,L,j,J,M> = (-l)^(L + 1)(L + 2) x 



J L 



JI+1 j + 1 



W + D \ lj + 1L U 1 j 



J L 



J L + 1 j - 1 



1 j — 1 L + ljl 1 j L + 2 



x 



Xn',L+2(r) 



( 2McS \ (V(r\ po ^ 4. (^ + 1 )( 2 ^ + 3) 2L + 3 d 



r dr 



Xn,h{ r ) dr ( 69 ) 



For L' = L — 2 they are: 



(n',L-2,j,J,M|(X-V) 2 -iv 2 |n,L,j,J,M) = (-1) yJ(L(L - 1) x 



J L 



J L — 1 j + 1 



J L j 

/2M cfr 



J L — 1 j — 1 
1 j L - 2 



x 



Xn',L-2(r) 



.V h 2 

And for the case V = L they are: 



(V(r) - < £ ) + 



L(2L-1) 2L-1 d 



r dr 



Xn,L(r)dr (70) 



(n', L, j, J, M| (A • V) 2 - 1 V 2 |n, L, j, J, M) = 



(L + l)(j + l) 
+W + 1) 



3 

J L j 
1 3+1 L+l 



+ (L + l)j 



R 



Xn'A r ) 



J L j 
1 j + 1 L - 1 

/2_Meff 

V ft 2 



+ A7 



J L j 
1 j - 1 L - 1 



J L j 
Ij-lHl 

2 



x 



(V(r)-< L ) Xn,L(r)dr 



For the case j = 0, J = L = L', and the matrix element reduces to: 
(ri, L, 0, L, M\ (X • V) 2 - ^V 2 |n, L, 0, L, M) = 0. 



(71) 



(72) 



which is the reason that the — |V 2 term was included as part of the hydrodynamic term. 
V(r) is the isotropic part of V (r, Xj and the basis states are assumed to be eigenstates of 
Ho defined by: 
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H » = -2^ V2 + V ' W (73) 

with eigenvalues E^ L . While not obvious by inspection, the above matrix elements can be 
shown to describe a Hermitian operator, as they should. 

Using this hamiltonian, the j — — > 1 spectrum of HCN in a He cluster can be recal- 
culated. Figure |IX| shows the spectrum calculated including the isotropic (Ceo) term in the 
potential and the hydrodynamic coupling, but without the anisotropic potential term (C 62 ) 
for a cluster of R = 3nm. The calculated spectrum of the impurity in the cluster looks 
remarkably like that (see Figure [IX]) for a linear molecule impurity in bulk liquid He, except 



the spectrum is 'inverted'. This can be understood as a consequence of averaging over the 
classical trajectory of the motion, which is in a plane perpendicular to L. For circular orbits, 
this leads to the prediction that the j '• = 1 state with rrij = is higher in energy than the 
states with rrij = ±1 by an amount equal to 2(— Cn)M e f[E^/5h 2 , where is the average 
kinetic energy for the center of mass motion. Using the thermal RMS kinetic energy, this 
gives a predicted average splitting of ~ 150 MHz. This naturally explains the shape of the 
predicted spectrum and the energy scale. 

Figure |Xj shows the calculated cluster spectrum for the HCN R(0) line including both 
the anisotropic potential term and the hydrodynamic coupling. Somewhat surprisingly, the 
spectrum is narrower than for the case of either coupling term by itself. The reason is that 
for the dominant AL = matrix elements, the two sets of matrix elements destructively 
interfere. This is despite the fact that the two coupling terms are complementary in that 
the hydrodynamic coupling is most important when the impurity is near the center of the 
cluster (and thus has the highest velocity), while the potential anisotropy is most important 
at the outer turning point of the motion. The quantitative effect of this interference on 
the predicted spectrum will be a strong function of the particular parameters used in the 
calculation, including the cluster size. 



VII. APPLICATION TO THE R(0) TRANSITION OF OCS 

We will now consider the application of the present model to the case of OCS as an 
impurity, partly because of the elegant studies of the spectrum of this molecule performed 
in Gottingen, and also because it has the narrowest lines presently observed for an impurity 
in 4 He clusters ||. In fact, the linewidths estimated above for the R(0) transition of HCN, 
while being well below that observed in the R(0) transition HCN in the first C-H overtone 
band, are similar to that of the R(0) line of the OCS z/ 3 fundamental. This suggests that 
the present model may be more accurate in this case. 

At the author's request, Joanna Howson and Jeremy Hutson have calculated long range 
coefficients for the He-OCS pair |7|. They find C 60 = 14.2 eVA 6 , C 62 = 2.44 eVA 6 , C n = 
14.4 eVA and C 73 = 13.9 eVA . Calculations were done with a number of basis sets, but 
the above results were done with a Dunning augmented-cc-pVTZ basis and are reported 
to have ~ 10% uncertainty. The calculations use SCF polarizabilities and employed the 
CADPAC6.3 program package. 

In the present paper, we will continue to use only the lowest order, Cg, terms, though the 
type of scaling done above indicates that the C71 term may make a substantial contribution 
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in the present case. Higgins and Klemperer jnj have calculated a series of points on the 



He-OCS surface at the MP4 level of theory, which they kindly made available to the author. 
From the position of the minimum as a function of angle between the He-OCS vector and 
the OCS molecular axis, we approximate the molecule as an ellipsoid with a = b = 3.4 A 
and c = 4.7 A, or f = 1.38. 

Figure |X| shows the R(0) transition of OCS observed in 4 He clusters, as reported by 
Grebenev et al. [14]]. It has a FWHM of 150 MHz, and is asymmetric with a broad shoulder 
going to lower frequencies. Also in Figure [IX| is plotted the R(0) transition for OCS, pre- 
dicted from the above parameters, and convoluted with a Lorentzian of 30 MHz FWHM to 
wash out most of the fine structure in the predicted spectrum. This is justified both because 
we expect vibrational and rotational relaxation lifetime to broaden the lines, and because 
the average over the cluster size distribution should at least partially wash out the high res- 
olution structure. The agreement with experiment is remarkable, especially considering the 
uncertainty in the parameters going into the model. The calculated linewidth is somewhat 
too small, but the shape of the peak is in excellent agreement with what was observed in 
the experiment. 

This level of agreement is highly suggestive that the interactions included in the present 
model are the principal source of line broadening in the OCS R(0) transition. However, 
it must be noted that the experiments have shown that the linewidth of the ro-vibrational 
lines increases dramatically with J, with the R(0) line being the narrowest. While the author 
has not calculated lineshapes for higher J transitions because of memory limitations on the 
personal computer on which these calculations have been performed, it is likely that the 
present model will not predict an increase in linewidth with J, but will more likely predict 
the opposite. Also the experiments show a significant difference in linewidth for the R(0) 
and P(l) lines, which should have the same lineshape (except for an inversion) in the present 
model, and in fact any model that does not include a vibrational dependence to the He- 
impurity interaction. We conclude that while the interactions considered in this paper likely 
play an important role in the lineshape of the R(0) transition, as yet unidentified interactions 
also play an important role, especially at higher J. 



VIII. HYDRODYNAMIC CONTRIBUTION TO THE EFFECTIVE MOMENT OF 

INERTIA 

Let us now consider the hydrodynamic effect on the effective moment of inertia for rota- 
tion of our symmetric top impurity. It has been established experimentally that molecules 
rotate in liquid He clusters with an effective rotational constant that is smaller than that 
of the same molecule in the gas phase. This implies an increase in the effective moment of 
inertia for rotation. For light impurities with large rotational constants, fractional change of 
the rotational constants is small. For example, in the case of H 2 0, the rotational constants 
in the cluster are almost identical with the gas phase ||. In contrast, heavy impurities with 
small gas phase rotational constants, such as (CHs^SiCCH have a rotational constant in 
the cluster of only 20% as large as in the gas phase |4[]. Some reduction in rotation con- 
stant is expected because of hydrodynamic effects, i.e., He must move out of the way of the 
rotating impurity, which contributes to the kinetic energy. Using the model of the impurity 
as an ellipsoid, we can directly calculate this increased moment of inertia from the size and 
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shape of the molecule. The velocity potential for such motion is given by Milne-Thomson. 
Only an integral over the surface of the ellipsoid is needed to calculate the hydrodynamic 
contribution to the rotational kinetic energy. This leads to the following result for the hy- 
drodynamic contribution to the moment of inertia for rotation about the axis perpendicular 
to the symmetry axis of the ellipsoid: 

A/ B = /'A(0 + M' 7l r c 2 om (74) 

/'= IjwV(l + f 2 ) (75) 

A( = - J (l-OVo-70) (?6) 

* } i+e 2(i-a-a+a(«o-7o) 1 1 

r com is the displacement of the impurity center of mass from the geometric center of the 
ellipsoid and I' is the value of Is for the ellipsoid if filled with liquid He and this provides 
the natural scale for the effect. The other symbols are as defined above. Figure shows a 



plot of A(£) as a function of £. For £ ~ 1, one can use the expansion A(£) ~ |(£ — l) 2 . . .. 

As an application of this expression, we will again consider OCS. Using the parameters 
that a = b = 3.4A and c = 4.7A estimated from the ab initio points []38|1 , and r com = O.lA, 
the calculated enhancement of the B moment of inertia is 9.33 uA , compared to a value of 
I B = 83.1 uA z ', or a predicted rotational constant of OCS in the cluster of 90% of the gas 
phase value. This is a small fraction of the observed reduction in B, where the observed B 
value in He clusters is 36% of the gas phase B value ||. 

Another model for the reduction of rotational constants is to assume a solvation shell of 
He atoms that rotate with the impurity. In the case of SF 6 , such a model seems to naturally 
explain the observed reduction in B ||, but does not naturally explain the observed value 
of the A rotational constant for (SF 6 ) 2 |§, which by such a model would be expected to be 
at least 1/2 of the B value for the SFg monomer in the cluster, while the observed value 
(based upon a spectral fit) is less than 1/4. In their recent paper on OCS in mixed 3 He 
and 4 He clusters, Grebenev et al. |J propose a model for the effective rotational constants 
of impurities in He with a position dependent normal fluid density around the impurity 
that rotates with it and thus contributes to the moment of inertia. However, they give 
no microscopic description of how this normal fraction could be defined, nor how it could 
be measured, thus this model appears to lack empirical content at least until it is more 
precisely defined. It is interesting to point out in this context that previous fully quantum 
microscopic calculations have found that He atoms in the first solvation shell about SF 6 
participate in exchange and are thus part of the superfluid ||39)| . These authors did find that 
at 0.625 K the superfluid fraction in a SF 6 He39, 0.67(7), was less than previously calculated 
for a He64 cluster, 0.9(1). It would be interesting for these calculations to be repeated for 
larger He clusters (at least enough to 'fill' the second solvation shell around the SF 6 ), and 
at a temperature near that found for experimental clusters in order to make better contact 
with the presently available spectroscopic experiments. 



IX. SUMMARY 

In this paper, we have considered the energetics and dynamics for the motion of impu- 
rities in spherical liquid 4 He clusters, which are treated as superfluid liquid droplets. It is 
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found that the effective potential for motion of the impurity is determined by the long range 
part of the impurity-He interaction. The isotropic part 'traps' the impurity near the center 
of the cluster. The anisotropic part results in a coupling of the rotation and center of mass 
motion of the impurity, which leads to inhomogeneous broadening in the predicted IR and 
microwave spectra. Calculated broadening values are on the order of 0.1-1 GHz, on the 
order of the narrowest lines observed in IR spectroscopy of impurities in liquid He. Hydro- 
dynamic effects are also found to play an important role by coupling molecular rotation to 
the velocity of center of mass motion. The inhomogeneous spectral effects of the distribu- 
tion of thermally excited ripplons are found to be much smaller and are likely negligible in 
determining linewidths, though they likely play an important role in determining the time 
required for the impurity and cluster to come into equilibrium. Direct hydrodynamic effects 
are found to make a minor contribution to the dramatic reductions in effective rotational 
constants often found for impurity species in He cluster spectroscopy. 

This work will surely not be the last word on this subject. The interactions considered 
here, while certainly present, do not appear to be sufficiently strong to explain the observed 
linewidths of many impurities, such as SFg and HCN. This suggest other, perhaps more 
subtle, physical effects are playing a role. New nonlinear spectroscopic experiments that 
use either pump-probe methods and/or spectral hole burning are clearly needed to sort out 
homogeneous and inhomogeneous effects in He cluster spectroscopy and to provide critical 
tests of simple models. Microscopic theory will also surely play an important role in helping 
to sort out the various physical effects that can play a role. 
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TABLES 



TABLE I. Thermodynamic averages for the motion of SFq and ripplons in a He cluster at 
T c = 0.38 K. R is the radius of the cluster, Q the partition function for the center of mass motion 
or the ripplons, < E > is the thermally averaged energy, expressed in units of kbT c , C v is the 
heat capacity in units of kb, S is the entropy in units of kb, <L 2 > is the mean squared value for 
the thermal angular momentum, < a 2 > is the mean squared value for the displacement of the 
impurity from the center of the cluster, and < n > is the mean number of quanta of excitation of 
the ripplon modes. 
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FIGURES 

SF Center of Mass Energy Levels in 
a R = 3 nm He Cluster 
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FIG. 1. Center of mass motion energy levels for SFq in a He cluster of radius 3 nm. All levels 
below 1 cm -1 are plotted. The ordinate is L, the total orbital angular momentum quantum number 
with states of increasing n, the number of nodes in the radial direction, increasing vertically in 
each column. 
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SF Center of Mass Energy Levels in 
a R = 3 nm He Cluster 
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FIG. 2. Same as figure 1, but with n as the ordinate and states of increasing L increasing 
vertically in each column. 
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SF Radial Distribution Function in R = 3 nm He Cluster 
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FIG. 3. Comparison of the classical and quantum Radial Distribution Functions, P(a), for SFq 
in a 3 nm He cluster. The Quantum distributions has been calculated from the set of all eigenstates 
with wavenumber less than 2 cm . Also plotted, is the isotropic potential, U(a), that confines the 
SFg in the cluster. 
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Potential for HCN in R = 3 nm He Cluster 
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FIG. 4. Comparison of isotropic (Po) an d anisotropic (P2) potentials for an HCN impurity in 
a Helium cluster of 3 nm radius. 
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SF in R = 3 nm He Cluster 
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FIG. 5. Calculated Spectral lineshape for the v% fundamental of SFg based upon the 
dipole-induced dipole spectral shift, which depends upon the position inside the cluster 
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HCN J = -> 1 Transition in R = 3 nm He Cluster 
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FIG. 6. Simulation of the HCN R(0) transition in a 3 nm radius Helium cluster, using the 
anisotropic potential discussed in the text, including coupling between the HCN rotation and its 
center of mass motion relative to the center of the He cluster. 
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R = 5 nm with HCN impurity 
Shift for J = -> 1 Transtion 
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FIG. 7. Calculated shift in HCN j = l,rrij = level as a result of the interaction with the 
surface ripplons, as a function of impurity displacement from the center of the cluster. Also plotted 
is the radial distribution function in arbitrary units to show region of thermal population. 



33 



Hydrodynamic Constants 
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FIG. 8. The hydrodynamic constants a\ and 71 and their difference as a function of £, which 
is the ratio of the size of the ellipsoid along the symmetry axis divided by the size perpendicular 
to this. 
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Calculated Lineshape of Linear Molecule in bulk liquid He 




FIG. 9. Calculated Spectrum for the R(0) and R(2) lines of a linear molecule in bulk liquid He. 
The ordinate is plotted in units of Ajik{,T c h , where An is defined in the text, T c the temperature 
of the liquid He, and the rest are standard symbols for physical constants. 
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FIG. 10. Calculated HCN R(0) line shape, including the isotropic term in the potential and 
the hydrodynamic coupling, for a He cluster of R = 3 nm. 



36 



HCN in R = 3 nm Cluster, both terms 
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FIG. 11. Calculated HCN R(0) line, including the isotropic and anisotropic terms in the po- 
tentials and the hydrodynamic coupling, for a He cluster of R = 3 nm. 
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R(0) Transition of OCS in 3 nm He Cluster 



H — I 

CO 
CD 




600 



200 

Frequency Shift / MHz 



600 



FIG. 12. Plot of R(0) line of the fundamental band of OCS in a 4 He cluster of mean size 
N = 2700 (R = 3.1 nm). Taken from the thesis of M. Hartmann. Also included is the calculated 
lineshape using anisotropic potential and hydrodynamic coupled, as decribed in the text, convoluted 
with a Lorentzian of 30 MHz FWHM. 
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Constant in Hydrodynamic Contribution to Moment of Inertia 
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FIG. 13. The hydrodynamic constant A(£) as a function of £. See Eq. [7^ for the definition of 
this quantity. 
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